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After a rapid increase in temperature across the deconfmement temperature Td to temperatures T S> Ik, pure 
gauge theories exhibit unstable long wavelength fluctuations in the approach to equilibrium. This phenomenon is 
analogous to spinodal decomposition observed in condensed matter physics, and also seen in models of disordered 
chiral condensate formation. At high temperature, the unstable modes occur only in the range < k < k c , where 
k c is on the order of the Debye screening mass mo- Equilibration always occurs via spinodal decomposition 
for SU(2). For SU(3) at T near Td, nucleation replaces spinodal decomposition as the dominant equilibration 
mechanism. Monte Carlo simulations of SU(2) and SU(3) lattice gauge theories exhibit the predicted phenomena. 
For SU(2), the observed value of k c is in reasonable agreement with a value predicted from previous lattice 
measurements of mo- 



1. Introduction 

Much of the experimental relevance of finite 
temperature QCD comes from non-equilibrium 
situations, as in the case of heavy ion collisions. 
Current lattice techniques are restricted to equi- 
librium simulations, and cannot address directly 
non-equilibrium time evolution. However, there 
are aspects of non-equilibrium behavior that are 
determined by the equilibrium effective potential. 
After a rapid change in temperature, induced by 
some external agency, it is the effective potential 
that determines the process of re-equilibration, in 
particular whether relaxation, nucleation or spin- 
odal decomposition is the dominant equilibration 
mechanism. We have shown that a large, rapid 
increase of temperature across the deconfmement 
temperature results in spinodal decomposition in 
pure SU{2) and SU(S) gauge theories. § In par- 
ticular, we have found the existence of exponen- 
tially growing long wavelength modes in the ap- 
proach to equilibrium, characteristic of spindodal 
decomposition. This behavior depends only 
on the features of the equilibrium effective action 
in an unstable region. 

We consider first the generic case of a pure 
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SU (N) gauge theory in which the temperature 
is raised rapidly from a temperature less than 
Td, the deconfmement temperature, to a temper- 
ature above Td- We refer to such a rapid increase 
in system temperature as a quench, a borrowing 
of nomenclature from condensed matter physics. 
The equilibrium order parameter for the decon- 
fmement transition is the Polyakov loop Lp. At 
temperatures T < Td, the Z(N) global symme- 
try associated with confinement is unbroken, and 
(Lp) = 0. When the temperature is rapidly in- 
creased to T > Td, (Lp) = is no longer the 
stable state of the system, and must evolve to a 
new equilibrium state with (Lp) =/= 0, reflecting 
the transition to the gluon plasma phase. 

2. Effective Potential 

The instability of (Lp) =0 at high temper- 
atures is seen directly from the one-loop finite 
temperature effective potential for the Polyakov 
loop, as derived by Gross, Pisarski and Yaffe 
and Weiss. @]g For simplicity, consider the case 
of SU(2). It is convenient to parametrize the 
Polyakov loop as 

L F = cos[7r(l-V)/2]. (1) 

where — 1 < ip < 1. The effective potential at 
one loop for gauge bosons in a constant Polyakov 
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loop background is 
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The one loop result dominates the effective po- 
tential for T » Aqcd due to asymptotic free- 
dom. The first term is the standard black body 
result, obtained when ip = 1. The use of the vari- 
able ip makes the Z(2) symmetry of the potential 
under ip — > —-0 manifest. Note that the equilib- 
rium value of ip is ±1, corresponding to Lp — ±1; 
ip = 0, corresponding to Lp = 0, is a maximum 
ofV(ip). 

Our picture of the quenching process is that 
the system is initially in a state where ip is equal 
to zero at some temperature below Td- When 
the system is quickly raised to a new temperature 
T > Td, the system is still in the state with ip = 0. 
However, the system is unstable, and must even- 
tually find its way to either ip = +1 or = — 1- 
Because we quench into a region of the phase dia- 
gram where V"(ip) < 0, the system will decay to 
the equilibrium state via spinodal decomposition. 

Because pure SU (2) gauge theory has a second- 
order deconfining transition, spinodal decomposi- 
tion will occur after quenching to any tempera- 
ture T > Td- The situation is more subtle for 
SU(3), where the deconfinement phase transition 
is first order. This implies the existence of a 
metastable confined phase for some range of tem- 
peratures above T^, which in turn implies that 
nuclcation is the dominant equilibration mecha- 
nism just above Td- Perturbation theory gives no 
hint of this behavior, and the one-loop effective 
potential is unstable at Lp = 0. For tempera- 
tures sufficiently large that the one-loop effective 
potential for L is reliable, spinodal decomposition 
will occur. 

3. Langevin Model 

In order to study the dynamics of this transi- 
tion, we use the effective action of Bhattacharya 
combined with Langevin dynamics. Q For SU(2), 
the effective action takes the form 



S eff [P] 
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The equilibrium distribution will be 
exp [— S e ff [ip]]. We postulate Langevin dynamics 
of the form 



dip(x, t) 



r){x,t) 



(4) 



dr 8ip(x 1 T / 

where the white noise r/ is normalized to 

(77(0;, T)r](x',T')) = 2TS 3 (x - x')5{t - t'). (5) 

The late-time relaxational behavior of ip is 
controlled by the Debye screening mass mo, 
given by m 2 D = 2g 2 T 2 /3. Near equilibrium, 
any effects of initial conditions decay exponen- 



tially as exp 
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early times, the initial conditions contribute to 
ip(k, T)ip(—k, t) ) a term of the form 



V;(fe,0)^(-fc,0)exp 
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where k 2 = g 2 T 2 /3 = m 2 D /2. Modes with k < k c 
are initially not damped but grow exponentially, 
with the k = mode growing the fastest. This is 
a characteristic feature of spinodal decomposition 
with a non-conserved order parameter. 

4. Monte Carlo Results 

We have verified the existence of spinodal de- 
composition with simulations of rapidly quenched 
pure SU (2) gauge theory. Lattices of size 32 3 x 4 
and 64 3 x 4 were equilibrated at /3 = 2.0, below 
the deconfinement transition at (3d = 2.2986 ± 
0.0006.0 The coupling constant was increased 
instantaneously to j3 = 3.0, and the approach to 
equilibrium monitored via the Polyakov loop and 
other observables. The heat bath algorithm was 
used both for equilibration at low temperature 
before the quench and for subsequent dynamical 
evolution after the quench. While this time evo- 
lution is not the true time evolution of the non- 
equilibrium quantum field theory, features such 
as spinodal decomposition which depend only on 
the equilibrium effective action will occur with 
any local updating algorithm which converges to 
the equilibrium distribution. The abrupt change 
in (3 is a potential cause of concern with this pro- 
cedure, since the lattice spacing, and hence the 
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physical volume, changes in all directions when /3 
is changed. However, the large spatial sizes used 
should mitigate this effect. 
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Figure 1. S(k,r) versus Monte Carlo time. 

Figure 1 shows the Fourier transform of the 
connected Polyakov two point function S(k,r) for 
low values of the wave number as a function of 
Monte Carlo time for the same simulation. Note 
the early exponential rise in these modes, fol- 
lowed by a sharp disappearance as the Polyakov 
loop reaches its equilibrium value, characteristic 
of spinodal decomposition. Only the low momen- 
tum modes exhibit this growth; above k c no such 
growth occurs. Although the general behavior 
of S(k,r) is the same for each run, many de- 
tails are run dependent. In this particular run, 
the k/T = 0.68 mode achieves a larger amplitude 
than k/T = 0.56, which is atypical. In some runs, 
there is clear evidence for mode-mode coupling, 
reflecting the nonlinearity of the system. For each 
run, we have estimated the rate of growth of each 
low-momentum mode by fitting log(S(k,T)) to a 
straight line in r for early times. We can extract 
k c as the value where the growth rate is zero. 
From equation (6), the growth rate of each line is 
proportional to k 2 — k 2 for the linearized theory, 
but this may not represent the true time evolution 
of the simulation. In any case, the growth rates 
measured in each run are highly dependent on ini- 
tial conditions. In figure 2, we plot these growth 
rates versus k 2 /T 2 for the same run used in figure 



1. The error bars are naive estimates of the error 
for each growth rate for this particular run. The x 
intercept provides an estimate of k 2 for each run. 
Using multiple 64 3 x 4 runs at /3 = 3.0, we have 
estimated k c /T to be 1.14 ± 0.02. The principle 
errors in this estimate come from the sensitivity of 
individual modes to initial conditions and the dis- 
crete character of k on the lattice. We can com- 
pare this with lattice measurements of the Debye 
screening length, assuming the one-loop relation 



fee _ m D (T) 
T ~ ^2T 
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holds in general. Using the results of Heller et 
a/.§ for m D (T), we obtain k c /T = 1.35(5). We 
consider this to be reasonable agreement, given 
the many uncertainties involved. 
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Figure 2. Growth rates versus k 2 /T 2 . 

It is natural to ask what happens when rapid 
cooling takes place, as might occur in the late 
stages of the expansion of a quark-gluon plasma, 
or in the early universe. At low temperatures, 
we expect that V{tp) has a single minimum and 
V"(i/j) > everywhere. Thus theory predicts the 
absence of unstable modes, and that relaxational 
processes should dominate the approach to equi- 
librium. Simulations of such a cooling process, in 
which 32 3 x 4 lattice configurations equilibrated 
at (3 = 3.0 are suddenly cooled to (3 — 2.0, shows 
no sign of spinodal decomposition. 
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We have begun studies of the more physically 
relevant case of SU(3). Since nucleation domi- 
nates at temperatures just above Td, it is nec- 
essary to go to higher temperatures to observe 
spinodal decomposition. Working with 32 3 x 4 
lattice configurations equilibrated at j3 = 5.5, we 
performed quenches to various values of (3 > (3d ~ 
5.69. For values of (3 close to (3d, we believe we 
see evidence for metastability. We have observed 
exponentially growing low wavelength modes for 
quenches to (3 > 5.80. In figure 3, we show the 
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Figure 3. S(k, r) versus Monte Carlo time. 

evolution of the three lowest non-zero modes as 
a function of simulation time for a quench from 
(3 = 5.5 to (3 — 5.92. Note the similarity to the 
behavior seen in SU(2). A systematic investiga- 
tion is underway, emphasizing k c as well as the 
limit of metastability. 

5. Conclusions 

Analytical and simulation results both indi- 
cate the relevance of spinodal decomposition in 
the equilibration of a gluon plasma after a rapid 
quench to high temperature. In the case of 
SU(N) gauge theories, we have shown that the 
physical parameter k c which controls domain 
growth can be determined from simulations. In 
the case of SU(3), determination of k c as a func- 
tion of (3 will allow us to define the limit of 
metastability as the point k c ((3 m ) = 0. Note that 



small values of k c require correspondingly large 
lattice sizes to see non-zero exponentially grow- 
ing modes. 

Just as the formation of a DCC|| may lead 
to enhanced production of low-momentum pions, 
spinodal decomposition could lead to enhanced 
production of low-momentum gluons in the early 
stages of plasma formation. The characteristic 
scale for such a phenomenon would be gT. 

Full QCD has quarks as well as gluons, and 
the phase structure is different. It would be 
very interesting to study the behavior of both 
the Polyakov loop and the chiral condensate af- 
ter abrupt changes in temperature. The Polyakov 
loop can be studied fairly easily if large equi- 
librated unqucnched lattice field configurations 
are available. However, the chiral condensate re- 
quires a reliable estimator for the local condensate 
as a function of space and simulation time. 
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